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Abstract 
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This semiannual progress report covers the period from September 1, 1986 
to February 28, 1986 under NASA grant NAG-1-441 entitled “Direct Solar-pumped 
Iodine Laser Amplifier". During this period the improvement on the collection 
system of the Tarmarack Solar-Simulator beam has been attemped. On the other 
hand the basic study of evaluating the solid state laser materials for the 
solar-pumping and also the work to construct a kinetic model algorithm for the 
fl ashl amp-pumped iodine lasers have been carried out. 

It was observed that the collector cone worked better than the lens 
assembly in order to collect the solar simulator beam and to focus it down to a 
strong power density. The study on the various laser materials and their 
lasing characteristics shows that the neodymium and chromium-co-doped 
gadolinium scandium gallium garnet (Nd:Cr:GSGG) may be a strong candiate for 
the high-power solar-pumped solid state laser crystal. On the other hand the 
improved kinetic modeling for the flash! amp-pumped iodine laser provides a 
good agreement between the theoretical model and the experimental data on the 
laser power output, and predicts the output parameters of a solar-pumped iodine 
1 aser. 
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I. Fesibility Study of the Solar-Pumped Dye Laser 


A. Introduction 

The parametric study of a dye laser amplifier pumped by a solar-simulator 
and a flash! amp, and the amplifier gain measurement at various pump-beam 
irradiances on a dye cell were reported in the previous semiannual report [Ref. 
1]. From the previous work it has been shown that a solar concentration of 
20,000 is required to reach the threshold of the rhodamine 6G dye laser. About 
maximum 5,000 solar constants was previously achieved at the focus of a 
cone-shaped collector of the front diameter of 9.6 cm with the Tarmarack solar 
simulator. Thus, a different approach to focus the solar-simulator beam has 
been carried out in this research period. 

An assembly of a couple of lenses was used instead of the conical 
collector to place the focus of the solar simulator, which is located in the 
shutter position, into the work region where the laser or laser amplifier can 
be assembled [refer to Fig. 1], Because of the reflection of the solar 
simulator beam on the lens surfaces, the solar simulator beam concentration, 
which was less than 1,000 solar constants, was not improved as compared to that 
obtained with the conical collector. Therefore, it was concluded that a solar 
simulator, which provides a higher power than the present Tarmarack 
solar-simulator does, or a larger solar furnace is necessary to provide an 
optical power density of at least 20,000 solar constants (=2,706 W/cm ) in 
order to attempt to make a solar-pumped dye laser. 

B. Experiment and Results 

Figure 1 shows the assembly of the several lenses to focus the solar 
simulator beam down to a small spot. The solar simulator beam at the focus of 
the lens assembly was reflected by the diffuse surface, which was made of a 
Taflon bar. The reflected beam was focused down by a f=12.5 cm lens onto the 
mochromator whose dial was set at 530 nm with a grating of 1180 grooves/mm and 
with a less than 5A linewidth. Then, the spectral intensity was detected by a 
photomultiplier tube. A silicon photodiode was used to pick up a reference 
signal and to trigger the oscilloscope. 

After the spectral irradiance of the solar simulator beam was measured. 


1 


Oscilloscope 
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Figure 1 Experimental setup to focus down the solar-simulator 

bam to a small cylindrical spot with the use of lenses. 
MC: monochrometer, PMT: photodiode. 
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the diffuse surface was turned by 180 and the standard lamp signal was 
observed in order to calibrate the solar simulator's intensity. Figures 2 and 
3 show the typically observed signals of the standard lamp and the 
solar-simulator beams, respectively, on a Nicolet oscilloscope. The observed 

_3 

spectral irradiance of the standard source was 1.45x10 solar constants 
2 

(=0.268 w/cm nm), and those of the solar simulator at various input currents 
are shown in Graph 1. The maximum intensity at the input current 600 A was 
about 777 solar constants. 

With an experimental setup as shown in Fig. 4 the intensity of the light 
at the focus (F3) of the three lens assembly was observed to be about 14% of 
that at FI. The light intensity at the focus (F2) of the first two lens 
assembly without the cylinderical lens was about 68% of that at FI. This 
results show that the more lens we use the more optical loss will be caused. 

In Fig. 5 the solar simulator beam collection was tried without the 
cylindrical lens. The observed intensity at the focus was about 512 solar 
constants which is lower than that obtained with the cylindrical lens. The 
focusing of the solar simulator beam became poor without the cylindrical lens. 

The spatial distribution of the solar-simulator beam along the focal line 
of the 9.6-cm conical collector was remeasured with the experimental setup 
shown in Fig. 6. The absolute intensity calibration was done in a two-step 
process. A sun-gun (UL model : SG-63A ) was placed instead of the standard lamp 
in Fig. 6 at a distance of 43.2 cm from the diffuser. The intensity of the 
sun-gun was calibrated with the standard lamp and a small number of neutral 
density filter. Then, this sun-gun was used to calibrate the intensity of the 
solar-simulator beam instead of using the standard lamp, because the standard 
lamp intensity was very weak to directly calibrate the solar simulator beam 
intensity and required a large number of ND filter. Possibly this may cause a 
large error. The measured absolute spectral irradiance of the sun-gun at 
*= 530 nm was 1.08 solar constants, and that of the solar simulator 5,204 solar 
constants at the maximum focus. 

In conclusion, the conical collector provides the higher concentration of 
the solar simulator beam than any other lens assembly does. The observed 
maximum concentration is still much lower than the required threshold for the 
dye laser, which is 20,000 solar constants. The study on the solar-pumped dye 
laser will be closed at this point until a large solar furnace or a high power 
sol ar-simul ator is available. 


3 


100 


1 ' 


LOWER BERM (mV) 


<3 

<3 to <3 in 

— • in- in cu 13 



<3 


in (3 in 

r\- in cu <3 

(A 1 **) utoe y3ddfl 


4 


Figure 2 A typical oscillogram of the standard lamp signal 

observed with the experimental setup shown in Fig. 
the monochrometer dial setting at 530nm. 
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Figure 3 A typical oscillogram of the solar simulator signal 

observed with the experimental setup shown in Fig. 1 at 
the monochromater dial setting at 530nm. The ND filter 
used was 5.0, and the solar-simulator input current 
600A. 








Figure 4 Experimental setup to test the light transmission 

through the lenses which are used to collect the solar 
simulator beam. 
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Figure 5 Two-lens assembly to collect the solar-simulator beam. 
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Figure 6 Experimental setup for the intensity measurement of the 
solar-simulator beam at the focus of a conical 
collector. The conical collector's size = 9.6 cm front 
I.D. x 5.6 cm rear I.D. x 8.9 cm width. 
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• II. Evaluation of the Solid State Laser Materials for the Solar-Pumping 

A. Introduction 

04 » 

Since the first observation of the optical maser action with the CaFg:Dy 
crystal by Kiss, Lewis and Duncan [Ref.l] in 1963, several kinds of the solid 
state materials, such as NdgOg-doped barium crown glass, Nd:YAG, Nd:Cr:YAG, 

Ruby and NdrCaWO^, have been used for the solar-pumped lasers. Among them, the 
Nd:YAG crystal has been known as an efficient laser materials, so far. Its 
highest laser output power observed is 100W with a solar pumping [Ref. 10]. 
However, the laser output power of at least an order of IkW cw is required for 
the space power transmission. In order to satisfy the demand, the improvement 
of the solar-pumped NdrYAG laser or the search of new laser material has to be 
carried out. This report shows the evaluation of the laser crystals such as 
Ruby, NdrYAG, Nd: Glass, NdrYLF, NdrCrrGSGG, Alexandrite and Emerald for the 
high power solar pumping. 

Section B will describe the historical background of the solar-pumped 
solid state lasers. In section C, the characteristics of the various solid 
state materials will be summarized. The temperature distribution in the 
crystals and the coolant flowrate will be determined in section D. Section E 
will show the experimental setup and the results on the spectral and spatial 
distribution of the Tarmarack solar-simulator beam at the focus of a conical 
collector. Finally, the laser system and the expected laser output will be 
discussed in section F. 

B. Historical Background of the Solar-Pumped Solid State Lasers. 

1963: Z.J. Kiss, H.R. Lewis, and R.C. Duncan in RCA Laboratories. [Ref.l] 

- Observed the first solar-pumped optical maser action with the 
CaF^-.Dy^ crystal at the liquid neon temperature of 27°K and at the 
wavelength of 2.3um. 

1964: G.R. Simpson in American Optical Company [Ref. 2] 

- Observed the laser action at the wavelength of 1.06^m and the 
temperature of 30°C with a 6.25 wt% Nd^Og-doped barium crown glass 
of 0.1 mm diameter and 30 mm long. 
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- pumping source = sun and carbon arc sun-simulator. 

- Flashlamp threshold = 1 Joule 

- Observed the laser oscillation at the input energies of 

2 

1 . 1 5 J ( = 144 w/cm = 1,062 solar constants) with sun-pumping and of 
2 

1.50J (= 188 w/cm = 1,386 solar constants) with the carbon arc 
sun-simulator pumping. 

1966: C.W. Reno in RCA Defense Electronics Products [Ref. 3] 

o 

- Measured the threshold input power of 100 W (= 35 w/cm = 259 
solar constants) for a solar-pumped Nd:YAG laser, and observed the 
laser output of 100 mW with a 25-W coupled power to the crystal. 

- Observed the same results for the Nd:Cr:YAG laser as those for the 
NdrYAG laser. The Cr-ions did not improve the laser action 
because of the slow Cr - Nd transfer time (=1 ms) compared to the 
fluorescence lifetime of the neodymium ions. 

1966: C.G. Young in American Optical Company [Ref. 4] 

- Observed that the threshold powers for the Ruby and Nd:CaW0 4 
cyrstals are larger than those for the Nd:YAG and Nd:Glass. 

- Observed a 1-W cw laser output from a Nd : YAG crystal of 3-mm 
diameter and 30-mm long and of 1 at % doping density with a 
sol ar-pumping. 

- Observed a 1 . 25-W pulse laser output from a sun-pumped Nd:Glass 
with a pulse width of 7 ms. 

1972: L. Huff in GTE Sylvania Inc. [Ref. 5] 

- Observed a 4.85 W multimode laser output from a Nd:YAG crystal 
with a 24-inch diameter solar collector. 

- Used a conductive cooling of the laser rod with a copper heat 
sink. 

- Observed a TEM qo output of 0.8-W. 

1976: J. Falk, L. Huff and J.D. Taynai in GTE Sylvania Inc. [Ref. 6] 

- Observed the mode-locked and internally frequency-doubled laser 
output from a Nd:YAG crystal 
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Multimode TEM CW TEM mode locked 
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- Used an end pumping on the crystal with a 24-inch diameter solar 
collector and with the use of a UV cut-off -fi Iter . 

1978: D. Radick, E. Reed, and C. Chadwick in GTE Sylvania [Ref. 7] 

- Developed a pulsed Nd:YAG laser for the space communication 
purpose with a 24-inch diameter collector. 

- Observed a TEM qo mode laser output of 400 mW. 

1981: V.M. Batenin, A.L. Golger, and I. I. Klimovski in USSR [Ref. 8]. 

- Theoretical study on the feasiblity of a solar-pumped color-center 
crystal laser. 

- The expected values for a waveguide-type crystal with a 

o 

solar-pumping of &10 solar constants: 

Gain <£10' 3 cm" 1 
efficiency £ 1-6% 

output powers 10-60 W from an 1-m active medium. 

1984: H. Arashi , Y. Oka, and N. Sasahara in Japan [Ref. 9]. 

- Observed an 18-W cw laser output from a Nd:YAG crystal of 4-mm 
diameter and 75-mm long and of 0.9 at % doping density with a 10-m 
aperture solar collector. 

2 

- Effective pump power = 1.1 kW( = 117 w/cm = 863 solar constants). 

- Conversion efficiency (pump power 1 aser power) = 1.64%. 

1986: M. Weksler, J. Shartz, and Weizmann in Israel [Ref. 10]. 

- Observed a 100-W laser output from a Nd:YAG crystal with a solar 
furnace. 

1986: L. Zapata at NASA Langley Research Center [Ref. 11] 

- Observed a 28-W laser output from a Nd:Glass fiber bundle with a 
solar-simulator pumping. 
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Laser Crystals for the High-Power Solar Pumping. 

(1 ) . Ideal Laser Crystals. 

The characteristics which an ideal laser crystal for the high-power 
solar pumping should have are listed as below: 

(a) Efficient pumping (or absorption) bands near the peak of the 
solar spectrum. 

(b) Long fluorescent lifetime to lower the threshold input power 
for a cw operation and to provide a high energy storage. 

(c) High operation temperature at least higher than the room 
temperature to facilitate the rejection of heat. 

(d) Good thermal conductivity and high thermal resistance for a cw 
operation. 

(e) Mechanically strong. 

(2) . Characteristics of various laser crystals. 

Table 1 shows the list of the characteristics of the various 
crystals obtained from Refs. 12-35. The advantage and disadvantage 
of each crystal are discussed as follows: 

Ruby 

(i) Disadvantages 

- 3-level laser system 

: requires a high threshold input power 
: gives a poor quantum efficiency 

- Small gain coefficient 

(ii) Advantages 

- Broad absorption bands 

- Visible laser output 

- High thermal conductivity 

- Long fluorescence lifetime 

- The lowest doping concentration is possible. 
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Table 1. characteristics of laser crystals 
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CHARACTERISTICS OF LASER CRYSTALS (concluded) 



Nd: YAG (Neodymium-doped Yittrium Aluminum Garnet) 

(i) Advantages 

- 4-level laser system 

- Good thermal conductivity 

- Good Quantum efficiency 

- High gain coefficient 

- Low threshold input power required 

- Commercially well developed (600 Watts laser output) 

- Relatively good slope efficiency 

- 100 Watts cw laser operation with the solar-pumping has been 
achieved already 

Nd:Gl ass 

(i) Advantages 

- 4-level laser system 

- Cheap! easy to fabricate to a high doping density 

- Relatively low threshold input power required 

- Large absorption of the solar beam at high doping densities 

(ii) Disadvantages 

- Poor thermal conductivity 

: but still good for a short pulse laser operation at a low 
repetition rate 

- Broad fluorescence line width 

: increases the laser threshold 

: permits a good energy storagy and shorter light pulses 

Nd: YLF (Neodymium-doped Yittrium Lithium Fluoride) 

(i) Advantages 

- 4-level laser system 

- Long fluorescence lifetime (Twice of Nd:YAG) 

: requires a low threshold input power 
: provides a good energy storage 

- Low emission cross section 

: provides a good energy storage 

- Low thermal Tensing effect, (smaller than NdrYAG) 
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Nd : Cr : GSGG (Neodymium- and chromium-codoped Gadolinium Scandium 
Gallium Garnet) 

(i) Advantages 

- 4-level laser system 

- Broad absorption bands over the solar spectrum due to the Cr 3+ 

ions 

- Effective energy transfer from the Cr 3+ -ions to the Nd 3+ -ions 

- Low threshold input power required 

- High slope efficiency 

- A strong candidate for the high power laser 

Alexandrite 

(i) Advantages . 

- Tunability (700-800 nm) 

- Good thermal conductivity 

(ii) Disadvantages 

- Broad fluorescence linewidth 

: Require a high threshold input power 


Emerald 

(i) Advantage 

- Tunability (720-842 nm) 

(ii) Disadvantage 

- Short fluorescence time and broad fluorescence linewidth 

: Requires a high threshold input power 

(3) Absorption of the Solar Spectrum by the Laser Crystals 

The absorption spectrum of various laser crystals are shown in Fig. 1 
through 5, and compared with the spectral irradiance of the air-mass-zero solar 
spectrum. The amount of absorption of the solar beam per unit area by laser 
crystals has been calculated for the input solar spectrum of 1 solar constant 
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Figure 1 The absorption spectrum of the Nd:YAG crystal at 
doping density of 1 at% compared to the spectral 
irradiance of the air-mass-zero solar spectrum. 
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Figure 2 The absorption spectrum of the NdrYLF (ft-pol ari zed) 

crystal with a doping density of 2 at% compared to the 
spectral irradiance of the air-mass-zero solar spectrum. 
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Figure 3 The absorption spectrum of the Nd: YLF (<T-pol ari zed) 
crystal with a doping density of 2 atomic percent 
compared to the spectral irradiance of the air-mass-zero 
solar spectrum. 
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Figure 4 The absorption spectrum of a Nd:Glass compared to the 
air-mass-zero solar spectrum. The glass composition: 
Si0 2 = 66wt%, Nd 2 0 3 = 5wt%, Na 2 0 3 = 16wt%, BaO = 5wt%, 
A1 9 0 o = 2wt%, and Sb o 0 Q = lwt%. 
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WAVELENGTH, nm 

Figure 5 The absorption spectrum of the Nd:Cr:GSGG crystal at a 

doping density of Cr = 1 xlO^ions/cm^ and Nd = 2 x 10 
3 

ions/cm compared to the air-mass-zero solar spectrum. 



with an incident angle of 45 as shown in Fig. 9. The absorbed solar beam 
intensity, I abs , in the crystal at a given spectral range f rom ^ toX 2 is 

f*2 

( 1 ) 


obtained from 
I 

or 


abs 


. I (M (1 - e" NS>(A) * d )d^ 

An 0 


abs 


>2 

A, V * 1 (1 


- £ i (A). d N . 
e abs ' )*d* 


where I Q (A) 


Spectral irradiance of the incident light (W/cm ) 

3 

Absorbing molecules (or ions) per cm 

o 

Molecular absorption cross section (cm ) 

Pass length of the pump beam in the crystal (cm) 

Absorption cefficient (cm~^) = NO(A). 

Then the ratio of the absorbed solar beam intensity to the total incident solar 
beam becomes 

I n (A)(l - e~^* (Akl )‘dx 
>1 0 


N 

*(*) 

d 

^abs 


(A) 


R = 


1 ! o 


(?»)• dA 


The absorbed solar beam intensity in the central portion of the crystal rod 
within a 1 .0-mm diameter is calculated from 


I 


center 



I (A) e -N ^*)( d/2 ‘ 0.5mm) 

o ' 




The sample data and calculated results are listed in appendix A, and Graph 1 
shows the results on the absorbed solar beam intensity by the various crystals 
at the various doping densities according to Eq. (1). The amount of the solar 
beam absorption is the largest in the Nd:Cr:GSGG crystal. The absorption is 
larger in the Nd:YAG crystal than in the Nd:Glass at the same doping density. 
However, a Nd:glass of a high doping density can be easily fabricated compared 
to any other crystals, and the solar beam absorption can be increased with a 
high doping density. Nevertheless, the doping density can not be increased too 
high because of the colli si onal quenching effect which reduces the fluorescence 
lifetime at high doping densities. 

Graph 2 shows the absorbed solar beam intensity at the rod center of the 
Nd:YAG within a 1.0-mm diameter for various doping densities and for various 
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DOPING DENSITY, x 10^20 lons/cm^3 

Graph 2 The absorbed solar beam intensity within the 1-mm 
diameter center of the various sizes of the Nd:YAG 
crystal as a function of the Nd-ion doping deni sty. 



rod diameters. The amount of the solar beam absorption at the rod center 
decreases as the rod diameter increases. Thus, the rod size is limited by the 
threshold pump power requirement at the rod center. 

(4) Threshold Input Power and Slope Efficiency 

Table 2, which is cited from Ref. 31, shows the population inversion 

densities at the threshold and the threshold input-power densities of various 

laser crystals. The Nd:YAG has the lowest values among Ruby, Nd: Glass, 

Alexandrite and Emerald. However, it has been recently observed [Ref. 36] that 

the threshold of the Nd:Cr:GSGG is about one half of that of the Nd:YAG for a 

flash! amp pumping, and its slope efficiency is twice of the NdrYAG's [Ref. 30]. 

In Ref. 17, the threshold of the Nd:YLF was observed to be 51 OW and that 

of the Nd:YAG to be 600W with a tungustan lamp pumping. These values 

2 

correspond to 251 solar constants (=340W/cm ) and 295 solar constants 
2 

(=40W/cm ), respectively. This means that the Nd : YLF laser crystal can be 
easily pumped with a low power solar simulator beam. 

D. Crystal Temperature and Coolant Flowrate 

For a cylindrical laser rod and cooling system as shown in figure 6, the 

radial temperature distribution can be determined by solving an one-dimensional 

heat conduction equation [Ref. 37]: 

T(r) = T(r ) + (A n /4K)(r 2 - r 2 ) 
ooo 

2 

where A =P / r L:rate of heat generated per unit volume in the crystal rod. 

P a =I a bs ( so1ar constant) (rod surface area) (1- rod surface 
reflectance): rate of heat dissipated by the rod. 
k: Thermal conductivity of rod. 

In steady state, 

P, = (Heat removed from the crystal surface by the coolant). 

a 

Thus, 

p a . h Dl L(T R - T f ) 

where h : Surface heat transfer coefficient, 

D-j : Rod diameter. 

T r = T(r Q ) : Rod surface temperature. 

Tp : Mean fluid temperature. 




Table 2. Population inversion densities and 

thresholds of various laser crystals. 


i*i- 

Lasing 

Material 

8irTo‘n 3 ^ v 

'■^ c x^ 

N^ = No/ 2 

N^/tv 
p c= T(J _ 

[W/cm 3 ] 

Nth^v 

■n 

[J/cm 3 ] 

Ruby 

3-level 

X = 694.3 nm 
Av = 11 cm -1 
t 0 = 3.10* 3 sec 

10 19 

960 

4600 

Nd 3 *:glass 

4-level 
X = 1.06 
Av = 200 cm -1 
t 0 = 3.10* 4 sec 

9.10 15 

5.6 

2.7 

Nd 3+ :YAG 
4-level 
X = 1.064 
Av * 6 cm -1 
t 0 = 5.5.10* 4 sec 

2.10 15 

0.7 

0.6 

Alexandrite 

4-level 

XA* = 720 nm 
Av = 1000 cm* 1 
to = 3.2.10* 4 sec 

1.6 x10 17 

140 

70 

Emerald 

4-level 

XA* * 740 nm 
Av =* 1500 cm* 1 
tq = 6.5.1 0~ 5 sec 

4.1 x 10 16 

170 

17 


N (i $, = population inversion density at threshold in the meta- 
stable state for a 4-level system. 

= population inversion density at threshold for a 3-level 
system. 

t 0 = decay time associated with radiative laser transition. 
t| = refractive index. 

Av = the width of the gain linewidth at room temperature. 

t c = cavity lifetime — the time at which the energy is lost in 

the laser cavity. 

= speed of light. 

= lasing wavelength. 


c 

X 
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jacket for the calculation of the temperature 
distribution in the rod. 
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The surface heat transfer coefficient can be determined as follows: 

flo 

Yw 


(i) for a laminar flow (900 N Rg 2,000) 


h_ = 0.02 


N 0.45 N 0.5 N 0.05 w — \0.8 


W-P, "Re y 'Vr 


) ( r 


(ii) for a turbulant flow ( 1 2,000^ N Rg ^ 220,000) 

0*2 0 33 n « 

h t = 0.02 N Re N pr I*) 0 ' 53 


where N Re =4m*/^ji(D 2 + D-| ) : Reynolds number, 

Npr = Cp^/Kfy : Prandtl number, 

N Gr = ^°2 ’ VVw 9’Y (T r " T p)3^ 2 : Grashof number, 

K : Thermal conductivity of the coolant, 
w 

D 2 : Inside diameter of the coolant jacket, 

JA. : Viscosity of the coolant, 

Cp : Specific heat of the coolant 
m* : Mass flow rate 
f : Density of coolant, 

: Volumetric thermal expansion coefficient, 
g : Gravitational constant. 

The axial temperautre gradient along the laser crystal can be approximated as 

AT R =at f =l Y m *V 

The sample data and calculations are listed in Appendix B. Graph 3 shows 
the radial temperature distribution of the various cyrstals with a rod shape of 
3.2mm diameter x 75mm long for an input solar beam intensity of 4,000 solar 
constants and for a water flowrate of 0.5 gal on per minute. The rod 
temperature of the Nd:Cr:GSGG crystal is the highest among the four crystals 
because of its largest absorption of the solar spectrum. Since the thermal 
conductivity of the Nd:Cr:GSGG is higher than that of the Nd:YlF, the actual 
temperature of the GSGG crystal will be very low for the same laser output as 
the NdrYAG's or NdrYLF's. Graph 4 shows the radial temperature distribution of 
the Nd:YAG of an 1.0 atomic percent doping density for various coolant 
flowrates. There is not a significant temperature difference at the coolant 
flowrates above 0.5 galon/minute. Graph 5 shows the axial temperature 
distribution of the Nd:YAG crystal of an 1.0 atomic percent doping density for 
various coolant flowrates. The flowrates below 0.5 galon/minute are the 
laminar flows and those above 2.0 galon/minute are the turbulent flows. The 
temperature difference is not significant at the flowrates above 0.5 
galon/minute. 
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Graph 3 Radial temperature distribution of various laser 

crystals at a water flowrate of 1.5 gal /min. The laser 
rod size considered was 3.2 mm (diameter) x 75 mm 
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Graph 4 Radial temperature distribution of a NdrYAG crystal at 
various water flowrates. The conditions considered were 
the same as those in Graph 3. 
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JISTANCE FROM COOLANT ENTRANCE, cm 

Graph 5 Axial temperature distribution of a NdrYAG crystal at 
various water flowrates. The same conditions as shown 
in Graph 3 were considered. 



E. Spectral and Spatial Distribution of the Tarmack Solar-Simulator Beam 

The experimental setup used to measure the spatial distribution of the 

solar-simulator beam at the focus of a conical collector is shown in figure 7. 
The dimensions of the conical collector and the coordinates for the spatial 
distribution are shown in figure 8. Graphs 6, 7 and 8 show the three 

dimensional intensity distribution of the solar-simulator beam and the location 
of the laser crystal. 

The experimental setup used to measure the spectral distribution of the 

solar-simulator beam at the focus of the conical collector is basically same as 

the one shown in figure 7 except the use of a lens between the spectrometer and 

the diffuse surface instead of the fiber optics. A quartz lens was used for 
the UV transmission. Graph 9 shows the absolute intensity of the Tarmarack 
solar-simulator at the simulator input current of 200 A. The intensity at the 
focus of the conical collector is about 1,000 times higher than that of the 
air-mass-zero solar spectrum. At the wavelength of about 830nm there appear the 
line spectra of the xenon gas. 

Graph 10 shows the spectral intensities of the solar-simulator beam at the 
various solar-simulator input currents. The peak intensity at I = 600 A is 
about 6 times higher than that at I = 200 A. 

F. Laser System and Expected Laser Output 

Figure 9 shows the schematic diagram of the experimental setup for the 
laser system, and figure 10 shows the detailed diagram of the laser crystal 
holder and the cooling water jacket. The mirror has a reflectivity of 
99.5% at = 1.06 m and a curvature of 0.3 m. The reflectivity of the mirror 
M 2 , which has a curvature of 0.3 m, is changed to several values for the 
maximum laser output. 

The Tamarack solar-simulator provides 4-kW optical beam power [Ref. 38] at 
the electrical input power of 45 kW which is obtained with the input current of 
600 A and the DC voltage of 75 volts. Assuming that there is no loss on the 
conical reflector and using the same optical conversion efficiency as the one 
achieved by H. Arash, et al [Ref. 9], the expected maximum laser output power 
can be calculated for a Nd:YAG crystal as follows: 

P Nd*YAG = (°-°164)(4kW) = 65.6 W. 

The slope efficiency of the Nd:Cr:GSGG crystal can be assumed to be twice 
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Figure 7 The systematic diagram for the spatial and spectral 

measurements of the solar-simulator beam at the focus of 
a conical collector. 
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Figure 8 Dimensions of the conical collector used to focus down 
the solar-simulator beam, and the coordinates for the 
spatial distribution measurement of the beam at its 
focus. 
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DISTANCE ALONG COLLECTOR CONE AXIS, cm 

ph 6. Spatial Distribution of the So 1 ar*-S1 mu 1 ator Beam at 
the Focus of the Collector Cone 
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amph 8. Sp at 1 a 1 Distribution of the Sol ar— Slmul ator Beam at 
the Focus of the Collector Cone. 
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ine absolute intensity of the solar-simulator spectrum 
at the input current of 200 A compared with the 
air-mass-zero solar spectrum. 
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Graph 10 The spectral irradiances of the solar simulator beam 
measured with an OMA at various solar-simul ator input 
currents. 
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Figure 9 Setup for the solid state laser experiments with 
Tarmarack solar-simulator pumping. 
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Figure 10 Laser crystal holder and water-jacket with the mirrors 
attached to both ends. 



higher than that of the Nd:YAG crystal according to Refs. 30 and 36. Then the 
expected laser output with a Nd:Cr:GSGG crystal will be 

P M , r rcrr = 2 x 65.6 W = 131.2 W. 

Nd:Cr:GSGG 

G. Conclusion 

Nd:YAG, Nd:Cr;GSGG and Nd:YLF crystals are good candidates for the 
solar-pumped CW laser material. Among them, Nd:Cr:GSGG has the largest 
absorption of the solar spectrum, and possibly can be scaled up to a high power 
laser. The thermal problem is not significant for the crystals except the 
Nd:Cr:GSGG up to 4,000 solar constants pump power with a flowrate less than 2 
galon/minute. The detailed calculation and experimental measurements should be 
done during the next project period to determine the maximum pump power for 
each crystal before the thermal problem becomes significant. 
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Appendix A. Absorption of the solar beam by various laser crystals, 


Nd : YA6 


SOLAR SPECTRUM 

ABSORPTION SPECTRUM OF LASANT 
INCIDENT ANSLE (in degree) ■ 


RANSE < nm ) 

115.0 - 1000000.0 

300.0 - 910.0 


INTE6RATED SOLAR INTENSITY <mW/cm*2> - 135.3 


ROD 

PUMPIN6 

DOPING 

ABSORBED 

RATIO 

ABSORBED 

DIA. 

LEN6TH 

DENSITY 

SOLAR INTENSITY 

(Abs. S.I./ 

POWER AT CENTER 

( cm > 

( cm ) 

(ions/cm' > 3 > 

< mW/cm“2 ) 

Integ. S.I. XX) 

dia*1 mm , ( mW/cm''2 ) 

.10 

.14 

E.50E+19 

1 .363 

1 .01 

1 .363 



1 .30E+20 

2.607 

1.93 

2.607 



1.43E+20 

2.843 

2.10 

2.843 



1 . 90E+20 

3.668 

2.71 

3.668 



2.G0E+20 

4.819 

3.56 

4.819 



2.8BE+20 

5.225 

3.86 

5.225 



3.25E+20 

5.816 

4.30 

5.816 



3.90E+20 

6.753 

4.99 

6.753 



5.20E+20 

8.481 

6.27 

8.481 



G.50E+20 

10.047 

7.43 

10.047 

.20 

.28 

G.50E+19 

2.607. 

1.93 

1 .300 



1.30E+20 

4.819 

3. 56 

2.391 



1 . 43E+20 

5.225 

3.86 

2.589 



1 .90E+20 

6.613 

4.89 

3.259 



2 . G0E+20 

8.481 

6.27 

4.146 



2.8GE+20 

9.125 

6.74 

4.447 



3.25E+20 

10.047 

7.43 

4.874 



3.90E+20 

11.483 

8.49 

5.528 



5.20E+20 

14.050 

10.38 

6.664 



G.50E+20 

16.304 

12.05 

7.626 

.30 

.42 

6.S0E+19 

3.754 

2.77 

1 .244 



1 . 30E+20 

6.753 

4.99 

2.212 



1 . 43E+20 

7.291 

5.39 

2.382 



1 . 90E+20 

9.100 

6.73 

2.944 



2 . G0E+20 

11.483 

8.49 

3.662 



2.8GE+20 

12.292 

9.09 

3.899 



3.2SE+20 

13.441 

9.93 

4.231 



3.90E+20 

15.21 1 

1 1 .24 

4.730 



5.20E+20 

18.320 

13.54 

5.569 



G.50E+20 

21 .003 

15.52 

6.257 

.50 

.71 

6.50E+19 

5.816 

4 . 30 

1.147 



1 . 30E+20 

10.047 

7.43 

1 .934 



1 . 43E+20 

10.780 

7.97 

2.066 



1 .90E+20 

13.202 

9.76 

2.488 



2.G0E+20 

16.304 

12.05 

3.002 



2.86E+20 

17.338 

12.81 

3.167 



3.25E+20 

18.793 

13.89 

3.394 



3.90E+20 

21 .003 

15.52 

3.727 



5.20E+20 

24.811 

18.34 

4.270 



G.50E+20 

28.025 

20.71 

4.699 

.70 

.99 

6.50E+19 

7.G40 

5.65 

1 .065 
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1 .30E+20 
1 . 43E+20 
1 . 90E+20 
2.G0E+20 
2.86E+20 
3.2SE+20 
3.90E+20 
5.20E+20 
6.50E+20 

1.00 1.41 E.50E+19 

1 . 30E+20 
1.43E+20 
1 .90E+20 
2.60E+20 
2.8SE+20 
3.25E+20 
3.90E+20 
5.20E+20 
6.50E+20 

1.50 2.12 6.50E+19 

1 . 30E+20 
1 .43E+20 
1 . 90E+20 
2.60E+20 
2.8SE+20 
3.25E+20 
3.90E+20 
5.20E+20 
6.50E+20 

2.00 2.83 B.50E+19 

1 .30E+20 
1 .43E+20 
1 .90E+20 
2.G0E+20 
2.86E+20 
3.25E+20 
3.90E+20 
5.20E+20 
6.50E+20 

3.00 4.24 B.50E+19 

1 . 30E+20 
1 . 43E+20 
1 .90E+20 
2.G0E+20 
2.8GE+20 
3.25E+20 
3.90E+20 
5.20E+20 
G.50E+20 


12.812 
13.G87 
1 G .547 
20.148 
21.335 
22.995 
25.495 
29.740 
33.291 

10.047 

1G.304 

17.338 

20.678 

24.811 

26.157 

28.025 

30.809 

35.449 

39.211 

13.441 

21.003 

22.221 

26.106 

30.809 

32.315 

34.385 

37.420 

42.352 

46.228 

16.304 
24.81 1 
26.157 
30.404 
35.449 
37.042 
39.211 

42.352 
47.347 
51 .176 

21.003 

30.809 
32.315 
36.983 

42.352 
44.006 
46 . 228 
49.378 
54.225 
57.802 

*•****•••** 


9.47 

1 .727 

10.12 

1 .834 

12.23 

2.169 

14.89 

2.567 

15.77 

2.692 

17.00 

2.862 

18.84 

3.109 

21.98 

3.504 

24.58 

3.808 

7.43 

.966 

12.05 

1 .498 

12.81 

1 .580 

15.28 

1 .834 

18.34 

2.127 

19.33 

2.218 

20.71 

2.340 

22.77 

2.515 

26.20 

2.786 

28.98 

2.984 

9.93 

.841 

15.52 

1 .238 

16.42 

1 .297 

19.30 

1 .475 

22.77 

1 .675 

23.88 

1 .736 

25.41 

1.816 

27.66 

1 .926 

31 .30 

2.085 

34.17 

2.188 

12.05 

.749 

18.34 

1 .063 

19.33 

1.108 

22.47 

1.244 

26.20 

1 .391 

27.38 

1 .434 

28.98 

1 .490 

31 .30 

1 .563 

34.99 

1 .659 

37.82 

1 .71 1 

15.52 

.619 

22.77 

.837 

23.88 

.867 

27.33 

.955 

31.30 

1 .042 

32.53 

1 .065 

34.17 

1 .093 

36.50 

1.126 

40.08 

1.156 

42.72 

1.158 

******** 

* * * 
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Nd'.YAG 


RANGE (nm) 

SOLAR SPECTRUM 115.0- 1000000.0 
ABSORPTION SPECTRUM OF LASANT 300.0 - 900.0 
INCIDENT ANGLE (in degree) * 45 


***********•**«***•»**•*•****«*«*** 
INTEGRATED SOLAR INTENSITY (mW/cm A 2> - 135.3 


ROD 

PUMPING 

DOPING 

ABSORBED 

RATIO 

ABSORBED 

DIA. 

LENGTH 

DENSITY 

SOLAR INTENSITY 

(Abs. S.I./ 

POWER AT CENTER 

( cm ) 

( cm ) 

( ions/cm A 3 ) 

( mW/cm A 2 ) 

Integ. S.I. XX) 

dia*1 mm , ( mU/cm': 

.10 

.14 

1 . 30E+20 

2.594 

1.92 

2.594 



1 . 43E+20 

2.829 

2.09 

2.829 



1 .95E+20 

3.734 

2.76 

3.734 



2.60E+20 

4.793 

3.54 

4.793 



3.90E+20 

6.714 

4.96 

6.714 

.20 

NJ 

00 

1 . 30E+20 

4.793 

3.54 

2.378 



1 .43E+20 

5.197 

3.84 

2.574 



1 .95E+20 

6.714 

4.96 

3.308 



2.60E+20 

8.429 

6.23 

4.121 



3.90E+20 

1 1 .407 

8.43 

5.490 

.30 

.42 

1 . 30E+20 

6.714 

4.96 

2.199 



1.43E+20 

7.248 

5.36 

2.368 



1 . 95E+20 

9.224 

6.82 

2.981 



2.G0E+20 

11 .407 

8.43 

3.636 



3.90E+20 

15.099 

11.16 

4.693 

.50 

.71 

1 . 30E+20 

9.983 

7.38 

1.922 



1 . 43E+20 

10.710 

7.92 

2.052 



1 .95E+20 

13.348 

9.87 

2.510 



2.S0E+20 

16. 181 

11.96 

2.978 



3.90E+20 

20.826 

15.39 

3.692 

.70 

.99 

1 . 30E+20 

12.724 

9.40 

1.715 



1 .43E+20 

13.591 

10.05 

1.820 



1 .95E+20 

16.700 

12.34 

2.183 



2.B0E+20 

19.982 

14.77 

2.543 



3.90E+20 

25.259 

18.67 

3.076 

* * * 

* * * * 

******* 

********* 

******** 

***** 
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Nd ; YLF CPHI_POLARIZED) 

RANGE Cnm) 

'SOLAR SPECTRUM 115.0 - 1000000.0 

ABSORPTION SPECTRUM OP LASANT 300.0 - 900.0 

INCIDENT ANGLE Cin degree) « 45 


x x x 

x x x x 

X X X X X 

X X 

X X X X X X 

xxxxxxxxx 

X X 

XXX 

X 

INTEGRATED SOLAR INTENSITY 

CmW/cm*2) * 

135.3 




ROD 

PUMPING 

DOPING 


ABSORBED 

RATIO 


ABSORBED 

D I A . 

Length 

DENSITY 


SOLAR INTENSITY CAbs. S.I./ 


POWER AT CENTER 

Ccm) 

Com) 

Cions/cm" 

‘3 ) 

CmW/cm-2) 

Integ. S.I. 

)CX> 

dia- 

1mm, CmW/cm* 

.10 

.14 

1 . 30E+20 


1.647 

1.22 



1.647 



2.SQE+20 


3.165 

2.34 



3.165 

.20 

.28 

1.30E+20 


3.165 

2.34 



1.580 



2.60E+20 


5.863 

4.34 



2.921 

.30 

. 42 

1 . 30E+20 


4.568 

3.38 



1.518 



2 . 60E+20 


8.208 

6.07 



2.705 

.50 

.71 

1 . 30E+20 


7.080 

5.23 



1.403 



2. 60E+20 


12.053 

8.91 



2.338 

.70 

.99 

1 . 30E+20 


9.262 

6.85 



1.301 



2 . 60E+20 


15.095 

11.16 



2.043 

1.00 

1.41 

1 . 30E+20 


12.053 

8.91 



1.168 



2 . 60E+20 


18.662 

13.79 



1.698 

1.50 

2.12 

1 . 30E+20 


15.761 

11.65 



.988 



2.60E+20 


22.992 

16.99 



1.300 

2.00 

2,33 

1 . 3QE+20 


18.662 

13.79 



.849 



2.60E+20 


26.165 

19.34 



1.039 

3.00 

4.24 

1 . 30E+20 


22.992 

16.99 



.650 



2.60E+20 


30.733 

22.71 



.728 

x x x 

x x x x 

* * * * * 

X X 

X X X X X X 

xxxxxxxxx 

X X 

XXX 

X X 
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•Nd : YLF CS I GMA_POLAR I ZED) 


SOLAR SPECTRUM 

ABSORPTION SPECTRUM OF LASAHT 
INCIDENT ANGLE Cin degree) - 


RANGE Cnm) 

115.0 - 1000000.0 

300.0 - 900.0 


***************** *xm****n****mMmmx* 
INTEGRATED SOLAR INTENSITY fmW/cm*2) - 135.3 


ROD 

PUMPING 

DOPING 

ABSORBED 

RATIO 

ABSORBED 

DIA. 

LENGTH 

DENSITY 

SOLAR INTENSITY 

CAbs. S.I./ 

POWER AT CENTER 

(cm) 

Ccm) 

Cions/c»“3 ) 

CmW/cn»*2) 

Integ. S. I.)CX> - 

dia*lmm, CmW/cm*2) 

.10 

.14 

6.50E*19 

.618 

.46 

.618 



1.30E+20 

1.216 

.90 

1.216 



1 . 95E+20 

1.794 

1.33 

1.794 



2. 60E+20 

2.355 

1.74 

2.355 



2. 36E+20 

2.574 

1.90 

2.574 



3.25E+20 

2.897 

2.14 

2.897 



3. 90E*20 

3.423 

2.53 

3.423 



5 . 20E+20 

4.427 

3.27 

4.427 



6. 50E+20 

5.372 

3.97 

5.372 

.20 

.28 

6. 50E+19 

1.216 

.90 

.608 



1.30E+20 

2.355 

1.74 

1.176 



1.95E+20 

3.423 

2.53 

1.709 



2 . 60E*20 

4.427 

3.27 

2.207 



2 . 86E+20 

4.811 

3.56 

2.397 



3. 25E*20 

5.372 

3.97 

2.674 



3. 90E+20 

6.263 

4.63 

3.112 



5. 20E+20 

7.900 

5.84 

3.908 



6.50E*20 

9.371 

6.93 

4.611 

.30 

.42 

6. 50E+19 

1.794 

1.33 

.598 



1.30E+20 

3.423 

2.53 

1.139 



1. 95E+20 

4.906 

3.63 

1.628 



2.60E+20 

6.263 

4.63 

2.072 



2.86E+20 

6.773 

5.01 

2.238 



3.25E+20 

7.508 

5.55 

2.474 



3.90E+20 

8.655 

6.40 

2.840 



5. 20E*20 

10.700 

7.91 

3.474 



6.50E+20 

12.473 

9.22 

4.000 

.50 

.71 

6.50E+19 

2.897 

2.14 

.579 



1.30E+20 

5.372 

3.97 

1.068 



1 . 95E+20 

7.508 

5.55 

1.483 



2.80E-20 

9.371 

6.93 

1.836 



2. 86E+20 

10.052 

7.43 

1.962 



3. 25E+20 

11.013 

8.14 

2.136 



3. 90E+20 

12.473 

9.22 

2.392 



5. 20E*20 

14.969 

11.06 

2.800 



6. 50E+20 

17.043 

12.60 

3.102 

.70 

.99 

6. 50E + 19 

3.932 

2.91 

.560 


• 

1 . 30E+20 

7.104 

5.25 

1.004 



1 . 95E+20 

. 9.716 

7. i8 

1.356 



2. 60E+20 

11.909 

8.80 

1.638 



2.86E»20 

12.692 

9.38 

1.734 



3. 25E *20 

13.784 

10.19 

1.863 



3 . 90E*20 

15.413 

11.39 

2.045 



5 . 20E*20 

13.135 

13.40 

2.314 



6 . 50E+20 


20.354 


15.04 


2.496 


1.C0 

, 1.41 

6.50E+19 

5.372 

3.97 

.534 

• 


1 . 30E+20 

9.371 

6.93 

.917 



1 . 95E+20 

12.473 

9.22 

1.195 



2 . 60E+20 

14.969 

11.06 

1.397 



2 . 86E+20 

15.842 

11.71 

1.462 



3. 25E+20 

17.043 

12.60 

1.547 



3. 90E+20 

18.811 

13.90 

1.659 



5.20E-20 

21.723 

16.06 

' 1.812 



6.50E+20 

24.079 

17.80 

1.909 

1.50 

2.12 

6. 50E+19 

7.508 

5.55 

.494 



1 . 30E+20 

12.473 

9.22 

.796 



1 . 95E + 20 

16.050 

11.86 

.984 



2.60E*20 

18.811 

13.90 

1.105 



2.36E+20 

19.760 

14.60 

1.141 



3.25E+20 

21.057 

15.56 

1.186 



3. 90E+20 

22 . 956 

16.97 

1.242 



5. 20E+20 

26.071 

19.27 

1.317 



6.50E+20 

28.598 

21.14 

1.369 

2.00 

2.83 

6.50E+19 

9.371 

6.93 

.459 



1 . 30E*20 

14.969 

11.06 

.698 



1 . 95E+20 

18.811 

13.90 

.829 



2 . 60E+20 

21.723 

16.06 

.905 



2 . 86E+20 

22.719 

16.79 

.926 



3.25E+20 

24.079 

17.80 

.953 



3. 90E+20 

26.071 

19.27 

.988 



5. 20E+20 

29.348 

21.69 

1.037 



6. 50E + 20 

32.012 

23.66 

1.077 

3.00 

4.24 

6.50E+19 

12.473 

9.22 

.398 



1.30E+20 

18.811 

13.90 

.552 



1. 95E+20 

22.956 

16.97 

.621 



2. 60E+20 

26.071 

19.27 

.658 



2.86E+20 

27.138 

20.06 

.669 



3. 25E+20 

28.598 

21.14 

.684 



3. 90E+20 

30.740 

22.72 

.705 



5. 20E+20 

34.274 

25.33 

.740 



6. 50E+20 

37.156 

27.46 

.770 

x * * 

m * x m 

X X X X X X 

xxxxxxxxxx 

X X X X X X X 

X I X I X X 



Nd sGLASS 


RANGE Cnm) 

SOLAR SPECTRUM 115.0 - 1000000.0 

ABSORPTION SPECTRUM OF LASANT 300.0 - 900.0 

INCIDENT ANGLE Cin degree) - 45 

GLASS COMPOSITION t S102 - 66 wt X, Nd203 » 5 wt X 

Na20 * 16 wt X, BaO > 5 wt X 

A1203- 2 wt I, Sb203 - 1 wt X 

*X***X**»*M9***MM*XX***MX*******X*X 
INTEGRATED SOLAR INTENSITY CmW/cm-2) - 135.3 


ROD 

DIAMETER 

(cm) 


.10 


.20 


.30 


.50 


.70 


1.00 


PUMPING DOPING ABSORBED RATIO ABSORBED 

LENGTH DENSITY SOLAR INTENSITY CAbs. S. I. / POWER AT CENTER 


Ccm) 

Cwt X Nd203) 

CmW/cm*2) 

.14 

1.00 

1.300 


2. 00 

2.534 


3.00 

3.703 


5.00 

5.891 


7.00 

7.877 


9.00 

9.692 


10.00 

10.542 

.28 

1.00 

2.534 


2.00 

4.826 


3.00 

6.907 


5.00 

10.542 


7.00 

13.612 


9.00 

16.244 


10.00 

17.426 

.42 

1.00 

3.708 


2.00 

6.907 


3.00 

9.692 


5.00 

14.307 


7.00 

17.988 


9.00 

21.012 


10.00 

22.337 

.71 

1.00 

5.891 


2.00 

10.542 


3.00 

14.307 


5.00 

20.065 


7.00 

24.325 


9.00 

27.669 


10.00 

29.100 

.99 

1.00 

7.877 


2.00 

13.612 


3.00 

17.988 


5.00 

24.325 


7.00 

28.825 


9.00 

32.285 


10.00 

33.751 

1.41 

1.00 

10.542 


2.00 

17.426 


3.00 

22.337 


5.00 

29.100 


7. 00 

33.751 


Integ. S.!.)(X) dia=lmm, CmW/cm*2) 


• 961 

1.300 

1-87 

2.534 

2.74 

3.703 

4.35 

5.891 

5.82 

7.877 

7.16 

9.692 

7.79 

10.542 

1.873 

1.267 

3.57 

2.409 

5.10 

3.440 

7.79 

5.217 

10.06 

6.680 

12.01 

7.890 

12.88 

8.416 

2.741 

1.235 

5.10 

2.292 

7.16 

3.199 

10.57 

4.651 

13.30 

5.735 

15.53 

6*552 

16.51 

6.884 

4.354 

1.174 

7.79 

2.081 

10.57 

2.785 

14.83 

3.765 

17.98 

4.376 

20.45 

4.768 

21.51 

4.911 

5.822 

1.118 

10.06 

1.898 

13.30 

2.447 

17.98 

3.119 

21.30 

3.478 

23.86 

3.682 

24.95 

3.751 

7.792 

1.040 

12.88 

1.665 

16.51 

2.047 

21.51 

2.443 

24.95 

2.620 


9.00 

10.00 


37.268 

38.744 


27.54 

28.64 


2.710 

2.738 


1.50 2.12 1.00 14. 

2.00 22 . 

3.00 27. 

5.00 34. 

7.00 39. 

9.00 42. 

10.00 44. 

2.00 2.83 1.00 17. 

2.00 26. 

3.00 31. 

5.00 38. 

7.00 43. 

9.00 46. 

10.00 48. 


XXXXXXXXXXXXXXXXXX 


307 

10.574 

.927 

337 

16.51 

1.364 

669 

20.45 

1.581 

715 

25.66 

1.763 

427 

29.14 

1.830 

927 

31.73 

1.860 

380 

32.80 

1.867 

426 

12.880 

.832 

088 

19.28 

1.143 

608 

23.36 

1.272 

744 

28.64 

1.366 

430 

32.10 

1.396 

864 

34.64 

1.403 

278 

35.68 

1.402 

i X X X 

xxxxxxxxxx 

X X X X 
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Nd : Cr : GSGG 


RANGE ( nm ) 

SOLAR SPECTRUM 115.0 - 1000000.0 

ABSORPTION SPECTRUM OF LASANT 300.0 - 920.0 

INCIDENT ANGLE (in degree) = 45 

DOPING DENSITY ( ions/cm' , 3 ) : Cr * 1.E+20 

Nd * 2.E+20 

*********************************** 

INTEGRATED SOLAR INTENSITY (mW/cm A 2) * 135.3 


ROD 

PUMPING 

DOPING 

ABSORBED 

RATIO 

ABSORBED 

DIA. 

LENGTH 

DENSITY 

SOLAR INTENSITY 

(Abs. S.I./ 

POWER AT CENTER 

( cm ) 

( cm ) 

( ions/cm'S ) 

( mW/cm , '2 ) 

Integ. S.I. XX) 

dia^l mm , ( mW/cm^C ) 

.30 

.42 

1 . 00E+00 

34.282 

25.34 

10.774 


*#****■*■»■** + •*•»*#**■#■»•***#*•*■**■*#******* 


Nd:Cr:GSGG 


SOLAR SPECTRUM 

ABSORPTION SPECTRUM OF LASANT 
INCIDENT ANGLE (in degree) * 


RANGE (nm > 

1 15.0 - 1000000.0 

300.0 ~ 900.0 


DOPING 

i DENSITY 

( ions/cm''3 ) : 

Cr = 

1 . E+20 






Nd = 

2.E+20 



* * * 

* * * * 

******* 

****** 

********* 

* * 

* * * * 

INTEGRATED SOLAR INTENSITY 

< mW/crTC ) * 

135.3 



ROD 

PUMPING 

DOPING ‘ 

ABSORBED 

RATIO 


ABSORBED 

DIA. 

LENGTH 

DENSITY 

SOLAR INTENSITY (Abs. S.I./ 


POWER AT CENTER 

( cm ) 

< cm ) 

< ions/cm'S ) 

( mW/cn A 2 ) 

Integ. S.I. 

>(%> 

dia=1 mm , ( mW/cm”2 ) 

.30 

.42 

1 .00E+00 

34 . 23G 

25.30 


10.758 

* * * 

* * * * 

******* 

****** 

********* 

* * 

***** 
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Appendix B. Temperature distribution in various laser crystals. 


CRYSTAL ****** 

CRYSTAL - 

ROD DIAMETER CcmD * 

ROD LENGTH CcmD = 

EFFECTIVE PUMPING LENGTH CcmD = 
THERMAL CONDUCTIVITY CW/cm/degree CD 


PUMPING SOURCE ****** 

TOTAL ABSORBED SOLAR INTENSITY CW/cm *2D * .006753 

SOLAR CONSTANT - 4000 

SURFACE ABSORPTION FACTOR - -8 

COOLING SYSTEM ****** 

NAME OF COOLANT - WATER 

THERMAL CONDUCTIVITY CW/cm/degree CD - .0569 

VOLUMETRIC THERMAL EXPANSION COEFF. Ci/degree 6.43E-5 

SPECIFIC HEAT CCal /g/degree CD = i 

CJ/g/degree CD * 4.186 

DENSITY C g/cm^D = 1 

VISCOSITY C g/cm/aec) = -01 

DIAMETER OF WATER JACKET CcmD = 1.308 

ENTERING TEMPERATURE OF COOLANT Cdegree CD = 20 

COOLANT MASS FLOW-RATE CGal /mi nD= .5 C- 31.6 cra*3/secD 


NdsYAG 

.32 

7.5 

7 

.14 
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PATE OF TOTAL HEAT DESSIPATION BY THE ROD (W/cm*2) « 152.07037576 

RATE OF HEAT GENERATION PER UNIT VOLUME IN THE CRYSTAL CW/cm^D - 270.12 

Prandt 1 Number * .735676625659 

Granshof Number without Temperature Term * 

Reynold Number » 1480.49291043 

COOLANT MASS FLOW-RATE CGal /mi nD - .3 C« 

LAMINAR FLOW ***** 

HEAT TRANSFER COEFFICIENT - 2.61334646526 

ROD TEMPERATURE AT THE COOLANT ENTRANCE Cdegree CD - 27.4775901055 

AXIAL TEMPERATURE GRAD I ANT OF ROD t COOLANT - 1.91908751961 


607.726091598 
18.9 cm~3/secD 


kxxxxkkxxkxxxkkkxkxxxkkxkxxkxk 

RADIAL TEMPERATURE DISTRIBUTION OF CRYSTAL 

# DISTANCE FROM ENTRANCE CcmD 3.75 


RADIAL DISTANCE FROM 
ROD CENTER CmmD 
0.000 
.160 
.320 
.480 
.640 
.800 
.960 
1.120 
1.280 
1.440 
1.600 


TEMPERATURE 

Cdegree CD 

40.7854767224 

40.6619932939 

40.2915430081 

39.6741258653 

38.8097418653 

37.6983910081 

36.3400732939 

34.7347887224 

32.8825372939 

30.7833190081 

28.4371338653 


X X X X X X X 


XKXKXXKXXXXXKXXXXXXXXXX 


RATE OF TOTAL HEAT DESSIPATION BY THE ROD CW/cm Jfc 2D - 152.07037576 

RATE OF HEAT GENERATION PER UNIT VOLUME IN THE CRYSTAL CW/cm*2D - 270.12 


Prandtl Number - .735676625659 

Granshof Number without Temperature Term * 607.726091598 

Reynold Number * 12337.4409203 

COOLANT MASS FLOW-RATE CGa 1 /mi nD - 2.5 C« 157.8 cm-3/secD 


TURBURANT FLOW * * * * 

HEAT TRANSFER COEFFICIENT * 4.11570935606 

ROD TEMPERATURE AT THE COOLANT ENTRANCE Cdegree CD * 25.2505165284 

AXIAL TEMPERATURE GRAD I ANT OF ROD i COOLANT - .230290502353 


XXKXXKXKXXXXXXXKXXXXXKKXXXXXXK 

RADIAL TEMPERATURE DISTRIBUTION OF CRYSTAL 

# DISTANCE FROM ENTRANCE Ccm) 3.75 


RADIAL DISTANCE FROM 
ROD CENTER CmmD 
0.000 
.160 
.320 
.480 
.640 
.800 
.960 
1.120 
1.280 


TEMPERATURE 

Cdegree CD 

37.7140046367 

37.5905212081 

37.2200709224 

36.6026537796 

35.7382697796 

34.6269189224 

33.2686012081 

31.6633166367 

29.8110652081 
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k 1.440 27.711 8469224 

1.600 25.3656617796 

*********************** ******* 


RATE OF TOTAL HEAT DESSIPATION BY THE ROD CW/cm*2! - 152.07037576 

RATE OF HEAT GENERATION PER UNIT VOLUME IN THE CRYSTAL CW/cm*2! - 270.12 

Prandt 1 Number - .735676625659 

Granshof Number without Temperature Term - 607.726091598 

Reynold Number « 2467.48818405 

COOLANT MASS FLOW-RATE (Gal /mi n!« .5 C- 31.6 cn»*3/sec! 

BOUNDARY Between LAMINAR and TURBULENT flows * * * * 

HEAT TRANSFER COEFF. Cosing LAMINAR flow equation! - 3.28873535843 

ROD TEMPERATURE at COOLANT ENTRANCE Cdegree C! - 26.0073594604 

HEAT TRANSFER COEFF. Cusing TURBULENT flow equation! « 1.1357132553 

ROD TEMPERATURE at COOLENT ENTRANCE Cdegree C! - 39.0273380179 

AXIAL TEMPERATURE GRAD I ANT OF ROD • COOLANT « 1.15145251177 


*K**KXXX**XX*X**K*KX*KXXSK*X*K 

RADIAL TEMPERATURE DISTRIBUTION OF CRYSTAL 

# DISTANCE FROM ENTRANCE Ccm! 3.75 


RADIAL DISTANCE FROM 
ROD CENTER Cmm! 

0.000 

.160 

.320 

.480 

.640 

.800 

.960 

1.120 

1.280 

1.440 

1.600 


TEMPERATURE 

Cdegree C! 

40.024069113 

39.9005856845 

39.5301353987 

38.9127182559 

38.0483342559 

36.9369833987 

35.5786656845 

33.973381113 

32.1211296845 

30.0219113987 

27.6757262559 


KKXXXXXXXXXXKXXXXXXXXXXXXXXXXX 


XXXXKXXXXXXXXXXXXXftXKKXXXXXXXX 

TEMPERATURE DISTRIBUTION ALONG THE LENGTH OF THE CRYSTAL 
AT CYSTAL SURFACE 


DISTANCE FROM TEMPERATURE 

ENTRANCE Ccm! Cdegree C! 


0.000 

<VJ 

.750 

27.2151452512 

1.5C0 

27.3302905024 

2.250 

27.4454357535 

3. 000 

27.5605810047 

3.750 

27.6757262559 


* 


4.500 

27.7908715071 

5.250 

* 27.9060167582 

6. 000 

28.0211620094 

6.750 

28.1363072606 

7.500 

28.2514525118 

* * X X X X X 

********* 


* 
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CRYSTAL ****** 

CRYSTAL * 

ROD DIAMETER Ccm) - 
ROD LENGTH Ccm) « 

EFFECTIVE PUMPING LENGTH Ccm) - 
THERMAL CONDUCTIVITY CW/cm/degree C) 


Nd a YLF Cl at X NdD 
.32 
7.5 
7 

.06 


PUMPING SOURCE ****** 

TOTAL ABSORBED SOLAR INTENSITY CW/cm*2) « .003996 

SOLAR CONSTANT - 4000 

SURFACE ABSORPTION FACTOR - .8 

COOLING SYSTEM ****** 

NAME OF COOLANT - WATER 

THERMAL CONDUCTIVITY CW/cm/degree C) - .0569 

VOLUMETRIC THERMAL EXPANSION COEFF. Cl/degree 6.43E-5 

SPECIFIC HEAT CCal /g/degree C) - 1 

CJ/g/degree C3 * 4.186 

DENSITY C g/c«n*33 • 1 

VISCOSITY C g/cm/sec) * .01 

DIAMETER OF WATER JACKET Ccm) - 1.308 

ENTERING TEMPERATURE OF COOLANT (degree C) - 20 

COOLANT MASS FLOW-RATE CGal /min)* .5 C« 31.6 cm*3/aec) 
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RATE OF TOTAL HEAT DESSIPATION BY THE ROD CW/cm *2) « 89.9856888192 

RATE OF HEAT GENERATION PER UNIT VOLUME IN THE CRYSTAL (W/cm*2) - 159.84 

Prandtl Number - .735676625859 

Granshof Number without Temperature Term * 

Reynold Number * 1480.49291043 

COOLANT MASS FLOW-RATE CGal /mi n) - .3 C 

LAMINAR FLOW ***** 

HEAT TRANSFER COEFFICIENT « 2.61334646526 

ROD TEMPERATURE AT THE COOLANT ENTRANCE (degree C) - 24.5367139586 

AXIAL TEMPERATURE GRAD I ANT OF ROD 4 COOLANT « 1.13559510267 


607.726091598 
18.9 ca*3/5ec) 


*K**XXX*KXKK**K*XXXXK*XKXKX*K* 

RADIAL TEMPERATURE DISTRIBUTION OF CRYSTAL 

# DISTANCE FROM ENTRANCE(cm) 3.75 


RADIAL DISTANCE FROM 
ROD CENTER (mm) 

0.000 

.160 

.320 

.480 

.640 

.800 

.960 

1.120 

1.280 

1.440 

1.600 


TEMPERATURE 

(degree C) 

42.15411151 

41.98361551 

41.47212751 

40.61964751 

39.42617551 

37.89171151 

36.01625551 

33.79980751 

31.24236751 

28.34393551 

25.10451151 


XXXXXXXXXKXKXXXKXKKXXXXXKKKXXK 


RATE OF TOTAL HEAT DESSIPATION BY THE ROD (W/cm*2) - 89.9856688192 

RATE OF HEAT GENERATION PER UNIT VOLUME IN THE CRYSTAL (W/cm*2) - 159.84 

Prandtl Number « .735676625659 

Granshof Number without Temperature Term - 607.726091598 

Reynold Number - 12337.4409203 

COOLANT MASS FLOW-RATE (Gal /mi n)- 2.5 (« 157.8 cm A 3/sec) 


TURBURANT FLOW * * * * 

HEAT TRANSFER COEFFICIENT - 4.11570935606 

ROD TEMPERATURE AT THE COOLANT ENTRANCE (degree C) - 23.1069249293 

AXIAL TEMPERATURE GRADIANT OF ROD 4 COOLANT « .136271412321 

XXXXXXXXXXXXXXXXXXXXKXKXKXXXXK 

RADIAL TEMPERATURE DISTRIBUTION OF CRYSTAL 

8 DISTANCE FROM ENTRANCE (cm) 3.75 


PAD I AL DISTANCE FROM 
ROD CENTER (mm) 

0.000 

.160 

.320 

.480 

.640 

.800 

.960 

1.120 

1.280 


TEMPERATURE 

(degree C) 

40.2246606354 

40 . 0541646354 

39.5426766354 

38.6901966354 

37.4967246354 

35.9622606354 

34.0868046354 

31.8703566354 

29.3129166354 
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' 1.440 26. 4144S46354 

1.600 23.1750606354 

kxxxxxkxkkxxxxxxkxkkkxxxxkxxxx 


RATE OF TOTAL HEAT DESSIPATION BY THE ROD (W/cm-2) - 89.9856688192 

RATE OF HEAT GENERATION PER UNIT VOLUME IN THE CRYSTAL (W/cro*2) « 159.84 

Prandtl Number » .735676625659 

Granshof Number without Temperature Term * 607.726091598 

Reynold Number « 2467.48818405 

COOLANT MASS FLOW-RATE (Gal/min)« .5 (• 31.6 cm*3/sec) 

BOUNDARY Between LAMINAR and TURBULENT flows * ■ * * 

HEAT TRANSFER C0EFF. (using LAMINAR flow equation) - 3.28873535843 

ROD TEMPERATURE at COOLANT ENTRANCE (degree C) - 23.644713221 

HEAT TRANSFER COEFF. Cusing TURBULENT flow equation) • 1.1357132553 

ROD TEMPERATURE at COOLENT ENTRANCE (degree C) - 31.2591800266 

AXIAL TEMPERATURE GRAD I ANT OF ROD * COOLANT - .681357061605 

KKKXSKXXKXKXKXKKXKKSKXXXXXXXXK 

RADIAL TEMPERATURE DISTRIBUTION OF CRYSTAL 

» DISTANCE FROM ENTRANCE (cm) 3.75 


RADIAL DISTANCE FROM 

TEMPERATURE 

ROD CENTER Own) 

(degree C) 

0.000 

41.7902785308 

.160 

41.6197825308 

.320 

41.1082945308 

.480 

40.2558145308 

.640 

39.0623425308 

.800 

37.5278785308 

.960 

35.6524225308 

1.120 

33.4359745308 

1.280 

30.3785345308 

1.440 

27.9801025308 

1.600 

24.7406785308 

XXXXXXXXXXX 

XXXXXXXK 


XXXKKKXXKKXXXXXXXXXXXXXXKXXXXX 

TEMPERATURE DISTRIBUTION ALONG THE LENGTH OF THE CRYSTAL 
AT CYSTAL SURFACE 

DISTANCE FROM TEMPERATURE 

ENTRANCE (cm) (degree C) 


0.000 

24.4 

.750 

24.4681357062 

1.500 

24.5362714123 

2.250 

24.6044071185 

3.000 

24.6725423246 

3.750 

24.7406785308 

4.500 

24.808814237 

5.250 

24.8769499431 

6.000 

24.9450856493 

6.750 

25.0132213554 

7.500 

25.0813570616 


XXXXXXXXXKKXKKKXXKKKXKKKXXXXKXX 
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CRYSTAL ****** 

CRYSTAL * 

ROD DIAMETER (cm) * 

ROD LENGTH (cm) - 

EFFECTIVE PUMPING LENGTH (cm) - 

THERMAL CONDUCTIVITY CW/c m/degree C) 


Nd. GLASS (ABOUT 1.3 wt 1) 
.32 
7.5 
7 

.0135 


PUMPING SOURCE ****** 

TOTAL ABSORBED SOLAR INTENSITY (W/cm-2) - .0045 

SOLAR CONSTANT « 4000 

SURFACE ABSORPTION FACTOR * .8 


COOLING SYSTEM ****** 

NAME OF COOLANT « WATER 

THERMAL CONDUCTIVITY CW/cm/degree C) « .0569 

VOLUMETRIC THERMAL EXPANSION COEFF. Cl/degree 6.43E-5 

SPECIFIC HEAT CCa 1 /g/degree C) « 1 

(J/g/degree C) * 4.186 

DENSITY C g/cm*3) « 1 

VISCOSITY ( g/cm/sec) » .01 

DIAMETER OF WATER JACKET (cm) - 1.308 

ENTERING TEMPERATURE OF COOLANT (degree C) • 20 

COOLANT MASS FLOW-RATE (Go I /min)- .3 C- 18.9 cm-3/sec) 


RATE OF TOTAL HEAT DESSIPATION BY THE ROD CW/cm*2) - 101.335212634 

RATE OF HEAT GENERATION PER UNIT VOLUME IN THE CRYSTAL CW/cm*2) - 180 

Prandtl Number « .735676625659 

Granshof Number without Temperature Term « 

Reynold Number * 1480.49291043 

COOLANT MASS FLOW-RATE (Gal /min)- .3 (« 

LAMINAR FLOW ***** 

HEAT TRANSFER COEFFICIENT - 2.61334646526 
ROD TEMPERATURE AT THE COOLANT ENTRANCE (degree C) « 25.0800958497 

AXIAL TEMPERATURE GRAD I ANT OF ROD * COOLANT • 1.27882331382 

XXXXXXXKXXXXXKXXXXXXKXXXXXXXKX 

RADIAL TEMPERATURE DISTRIBUTION OF CRYSTAL 

# DISTANCE FROM ENTRANCE (cm) 3.75 


607.726091598 
18.9 cm*3/sec) 


RADIAL DISTANCE FROM 

TEMPERATURE 

ROD CENTER (mm) 

(degree C5 

0.000 

111 . 05284084 

.160 

110.199507507 

.320 

107.639507507 

.480 

103.37284084 

.640 

97.3995075066 

.800 

89.7195075066 

.960 

80.3328408399 

1.120 

69.2395075066 

1.280 

56.4395075066 

1.440 

41.9328408399 

1.600 

25.7195075066 

xxxxxxxxxxx 

XXXKXXXK 




RATE OF TOTAL HEAT DESSIPATION EY THE ROD CW/cm»2) - 101.335212634 

RATE OF HEAT GENERAT I OH PER UNIT VOLUME IN THE CRYSTAL CW/cm-2) - 180 

Prandtl Number - .735676625659 

Granshof Number without Temperature Term « 607.726091598 

Reynold Number * 12337.4409203 

COOLANT MASS FLOW-RATE CGal /min)- 2.5 (- 157.8 cm*3/5ec) 

TURBURANT FLOW x X m * 

HEAT TRANSFER COEFFICIENT - 4.11570935606 

ROD TEMPERATURE AT THE COOLANT ENTRANCE (degree C5 - 23.4987893348 

AXIAL TEMPERATURE GRAD I ANT OF ROD « COOLANT - .153458797659 

XXXXXXXXKXXXXKXKXXXKXXKXXXKXXX 

RADIAL TEMPERATURE DISTRIBUTION OF CRYSTAL 

« DISTANCE FROM ENTRANCE (cm) 3.75 


RADIAL DISTANCE FROM 

TEMPERATURE 

ROD CENTER 

Cdesree CD 

0.000 

108.908852067 

.160 

108.055518734 

.320 

105.495518734 

.480 

101.228852067 

.640 

95.2555187338 

.800 

87.5755187336 

.960 

78.1888520669 

1.120 

67.0955187336 

1.280 

54.2955187336 

1.440 

39.7888520669 

1.600 

23.5755187336 

* xxxxxxxxxx 

XXKKXKXX 
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RATE OF TOTAL HEAT DESS1PATION BY THE ROD (W/cm*2) * 101.335212634 

RATE OF HEAT GENERATION PER UNIT VOLUME IN THE CRYSTAL (W/cm*2) - 180 

Prandtl Number * .735676625659 

Granshof Number without Temperature Term *’ 607.726091598 

Reynold Number « 2467.48818405 

COOLANT MASS FLOW-RATE CGa 1 /mi n) * .5 (« 31.6 cm A 3/sec) 


BOUNDARY Between LAMINAR and TURBULENT flows * * * * 

HEAT TRANSFER COEFF. (using LAMINAR flow equation) « 3.28873535843 

ROD TEMPERATURE at COOLANT ENTRANCE (degree C) - 24.0812563181 


HEAT TRANSFER COEFF. (using TURBULENT flow equation) - 1.1357132553 

ROD TEMPERATURE at COOLENT ENTRANCE (degree C) - 32.6792567867 


AXIAL TEMPERATURE GRAD I ANT OF ROD t COOLANT - .767293988294 


XXXKXXXXXXKXXXXKXXXXXXXXXXXXXX 

RADIAL TEMPERATURE DISTRIBUTION OF CRYSTAL 

# DISTANCE FROM ENTRANCE (ex) 3.75 


RADIAL DISTANCE FROM 

TEMPERATURE 

ROD CENTER (mm) 

(degree C) 

0.000 

110.616980327 

.160 

109.763646994 

.320 

107.203646994 

.480 

102.936980327 

.640 

96.9636469941 

.800 

89.2836469941 

.960 

79.8969803275 

1.120 

68.8036469941 

1.280 

56.0036469941 

1.440 

41.4969803275 

1.600 

25.2836469941 

xxxxxxxxxxx 

xxxxxxxx 


xxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

TEMPERATURE DISTRIBUTION ALONG THE LENGTH OF THE CRYSTAL 
AT CYSTAL SURFACE 


DISTANCE FROM TEMPERATURE 

ENTRANCE (cm) (degree C) 


0.000 

.750 

1.500 

2.250 
3. 000 
3. 750 

4.500 

5.250 
6.000 
6.750 

7.500 


24.9 

24.9767293988 

25. 0534587977 

25.1301881965 

25.2069175953 

25.2836469941 

25.360376393 

25.4371057918 

25.5138351906 

25.5905645895 

25.6672939883 


xxxxxxxxxxxxxxxxxxxxxxxx 
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CRYSTAL ****** 

CRYSTAL - 

ROD DIAMETER Ccm) - 
ROD LENGTH Ccm> * 

EFFECTIVE PUMPING LENGTH Ccm) « 
THERMAL CONDUCTIVITY CW/cm/degrte C) 


Nd « Cr sGSGG CNd-2E20D 
.32 
7.5 
7 

.08 


PUMPING SOURCE ****** 

TOTAL ABSORBED SOLAR INTENSITY CW/cm*2) 

SOLAR CONSTANT - 

SURFACE ABSORPTION FACTOR - 


.034282 

4000 

.8 


COOLING SYSTEM ****** 

NAME OF COOLANT - WATER 

THERMAL CONDUCTIVITY CW/cm/degree C) - .0569 

VOLUMETRIC THERMAL EXPANSION COEFF. Cl/degree 6.43E-5 

SPECIFIC HEAT CCal /g/degree C> « 1 

CJ/g/degree C) - 4.186 

DENSITY < g/cm*3) - 1 

VISCOSITY C g/cm/sec) - .01 

DIAMETER OF WATER JACKET Cc») - 1.308 

ENTERING TEMPERATURE OF COOLANT Cdegree C> - 20 


COOLANT MASS FLOW-RATECGa 1 2.5 C- 157.8 cm*3/sec> 



% 

RATE OF TOTAL HEAT DESSIPATJON BY THE ROD (W/cm*2) « 771-994168783 

RATE OF HEAT GENERATION PER UNIT VOLUME IN THE CRYSTAL CW/cm*2D - 1371-28 

Prandt 1 Number * .735676625659 

Granshof Number without Temperature Term * 

Reynold Number - 1480.49291043 

COOLANT MASS FLOW-RATECGal /min) « .3 O 

LAMINAR FLOW ***** 

HEAT TRANSFER COEFFICIENT - 2.61334646526 

ROD TEMPERATURE AT THE COOLANT ENTRANCE (degree CD • 55.134397703 

AXIAL TEMPERATURE GRAD I ANT OF ROD * COOLANT • 9.74236018766 


607.726091598 
18.9 cm*3/secD 


KXXXXXXKXXXKXXXXXXKXKXXXXXMXKK 

RADIAL TEMPERATURE DISTRIBUTION OF CRYSTAL 

# DISTANCE FROM ENTRANCE (cmD 3.75 


RADIAL DISTANCE FROM 
ROD CENTER CmmD 
0.000 
.160 
.320 
.480 
.640 
.800 
.960 
1.120 
1.280 
1.440 
1.600 


TEMPERATURE 

(degree CD 

169.707977797 

168.610953797 

165.319881797 

159.834761797 

152.155593797 

142.282377797 

130.215113797 

115.953801797 

99.4984417968 

80.8490337968 

60.0055777968 


XXXXXXXXXXXXXXXXXXXXXXXXXXXXKX 


RATE OF TOTAL HEAT DESSIPATION BY THE ROD CW/cm < *2) - 771.994168783 

RATE OF HEAT GENERATION PER UNIT VOLUME IN THE CRYSTAL CW/cm^D - 1371.28 

Prandtl Number « .735676625659 

Granshof Number without Temperature Term » 607.726091598 

Reynold Number - 12337.4409203 

COOLANT MASS FLOW-RATE (Ge 1 /miiO- 2.5 (« 157.8 cm*3/secD 

TUR8URANT FLOW * * * * 

HEAT TRANSFER COEFFICIENT « 4.11570935606 

ROD TEMPERATURE AT THE COOLANT ENTRANCE (degree CD « 46.654554661 

AXIAL TEMPERATURE GRAD I ANT OF ROD * COOLANT « 1.16908322252 


XX XXXXKXXXXXXXKXXXKXXXKXXXXXXX 

RADIAL TEMPERATURE DISTRIBUTION OF CRYSTAL 

# DISTANCE FROM ENTRANCE (cmD 3.75 


RADIAL DISTANCE FROM 
ROD CENTER CmmD 
0.000 
.160 
.320 
.480 
.640 
.800 
.960 
1.120 
1.280 


TEMPERATURE 

(degree CD 

156.941496272 

155.844472272 

152.553400272 

147.068280272 

133.389112272 

129.515896272 

117.448632272 

1C3. 18732027 2 

86.7319602723 
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RATE OF TOTAL HEAT DESSIPATJON BY THE ROD CW/cm*2) - 771.994168783 

RATE OF HEAT GENERATION PER UNIT VOLUME IN THE CRYSTAL CW/cm-2) - 1371.28 

Prandtl Number * .735676625659 

Granshof Number without Temperature Term * 607.726091598 

Reynold Number * 2467.48818405 

COOLANT MASS FLOW-RATE (Gal /mi n) * .5 (- 31.6 cm*3/sec) 


BOUNDARY Between LAMINAR and TURBULENT flows * * * * 

HEAT TRANSFER COEFF. (using LAMINAR flow equation) - 3.28873535843 

ROD TEMPERATURE at COOLANT ENTRANCE (degree C) - 48.2263340796 


HEAT TRANSFER COEFF. (using TURBULENT flow equation) - 1.1357132553 

ROD TEMPERATURE at COOLENT ENTRANCE (degree C) - 116.593395813 


AXIAL TEMPERATURE GRAD I ANT OF ROD & COOLANT - 5.8454161126 


XKXXXXXXXXXXXXXXXXXXXXXXXKXXXX 

RADIAL TEMPERATURE DISTRIBUTION OF CRYSTAL 

• DISTANCE FROM ENTRANCE (cm) 3.75 


RADIAL DISTANCE FROM 
ROD CENTER CmrrO 
0.000 
.160 
.320 
.480 
.640 
.800 
.960 
1.120 
1.280 
1.440 
1.600 


TEMPERATURE 

(degree C) 

167.025108056 

165.928084056 

162.637012056 

157.151892056 

149.472724056 

139.599508056 

127.532244056 

113.270932056 

96.8155720563 

78.1661640563 

57.3227080563 


XXXXXXXXXXXXXXXXXXKXKX 


xxxxxxxx 


xxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

TEMPERATURE DISTRIBUTION ALONG THE LENGTH OF THE CRYSTAL 
AT CYSTAL SURFACE 


DISTANCE FROM 
ENTRANCE (cm) 


TEMPERATURE 
(degree C) 


0.000 

.750 

1.500 

2.250 
3. 000 

3.750 

4.500 

5.250 

6.000 

6.750 

7.500 


54.4 

54.9845416113 

55.5690832225 

56.1536248338 

56.738166445 

57.3227080563 

57.9072496676 

58.4917912788 

59.0763328901 

59.6603745013 

60.2454161126 


XKXXXXXXKXXXKXXXXXKXXXK 
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III. Kinetic Modeling of Solar-Pumped Iodine Laser 


This report relates the work performed during the period of Sept. 1, 1986 
through February 28, 1987. During this time, the work to construct a kinetic 
model algorithm which predicts the output parameters of a solar-pumped iodine 
laser (lasing time, time to threshold, and pulse energy) has progressed. To 
this date the kinetic model has improved such that there is good agreement 
between the theoretical model and the experimental data for the system, defined 
in the previous progress report as flash! amp-pumped iodine laser oscillator 
system (characterized by a cylindrical resonator cavity and pumping times of 
about one millisecond). The three experimental data sets given graphically in 
Fig. 1, 2, and 3 are for three iodides i-C^F^I, n-C^Fgl, and t-C^Fgl used as 
lasants in this experimental setup. The subsequent kinetic rates found for the 
different gases are given in Table I. These kinetic rates are important to 
establish scalability for space to space power transmission. By introducing the 
heat equation, a diffusion time constant, the pumping rates, and the 
photodissociation cross-section for the three iodides (May 1986 Progress 
Report) into the kinetic model (NASA TP-2241); the model has been modified to 
incorporate a cylindrical laser cavity. 

Of the three lasants the rate coefficients have been best identified for i 

- CgFyl; and as reported previously with in experimentally found bounds for 
these coefficients, a best fit is found for the data given by this experiment. 
To find the best fit to the experimental data by the theoretical model an 
initial I 2 density is introduced as a function of fill pressure (Oct. 1986 
Progress Report) for the three lasants. In Fig. 2 and 3 for the gases n - 
C4F9I and t - C4F9I, the initial I 2 density is retained and C^, C^, and Q 2 
(Table I) are kept constant in the theoretical model, thereby a best fit to the 
data using the information known about the two iodides is found by adjusting 
the rate coefficients. Therefore, the reaction rates are identified in the 
same system with similar initial conditionsgiving the comparative reaction 
rates which are tabulated in Table I. 

Of the three iodides the least is known about the rate coefficients for n 

- C^Fgl. In spite of this, very good agreement has been obtained between the 
theoretical predictions and the experimental data (Fig. 2) for this gas. The 
predicted lasing times are shorter than that given by experiment, but Fig. 7, 
8, and 9 show good agreement with the theoretical power output and the 
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experimental data for the power output, except very late in the pulse. The 
difference very late in the pulse implies a pumping mechanism which remains to 
be included in the model, wall effects are being investigated. 

The iodide t - C^Fgl has been identified as the preferable lasant 
material. As the reaction rates demonstrate (Table I), this lasant recycles 
easier since the recombination rate to the parent molecule is higher and the 
formation of the radical dymer is relatively small. In Fig. 3 the time to 
threshold predicted by the model is later than that of the experimental data 
for the gas t - C^Fgl. This implies that there is an omission in theoretical 
model using the rate coefficients found thus far. On ther hand, as shown in 
Fig. 10, 11 and 12, the power output as a function of time agrees with the 
experimental data as the pulse decays. This problem is being investigated and 
will be reported later. (Probably initial density is too large since the 
system maybe cleaner.) 

To identify the relative contribution of a specific reaction upon the gain 
in the medium, the derivative of the inversion density as a functionof pulse 
duration divided by the sum of the photodissociation rates is given in Fig. 13 
through 21 for the three gases and or different fill pressure. For instance, 
for the reaction 


*1 

I* + R — ^ RI 
d I/dt = -K 1 [I*] [R] 

Where the negative sign indicates a loss of inversion density, where as, a 

positive sign would have indicated a pumping reaction. The relative peak value 

for each rate is given in the key along with the reaction's sign (+ or -). If 

-5 

the contribution is bleow 10 , it is omitted from the graph. Also shown is 

the theoretical prediction for the power output normalized to one and given in 

2 

the key is the peak power (W/cm ). These graphs show explicitly that quenching 
by I^ turns off the laser - a reaction that increases at higher pressures. 
Furthermore, the graphs give the relative size of the contribution of each of 
the reactions to the inversion density. 

Upon examination of Fig. 1, 2 and 3, it can be seen that by adjusting the 
rate coefficients within the known bounds as reported previously in literature 
good agreement is found for the theoretical prediction of experimental data. 
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Therefore, a comparison of the important rate coefficients is accomplished 
(Tabel I) for these three iodides in the experiment design stated above. In 
addition, the theoretical model as found is used to indicate graphically the 
individual contribution by each of the chemical reactions to the inversion 
density of the lasant, showing the relative magnitudes and times the reaction 
contributes to the output power of the laser. 

It remains to simulate multiple pumping pulses using the same gas fill for 
different pressures and comparing the theoretical results to the experimental 
data. This would further characterize the Ig diffusion and the wall reactions 
in addition to be her understanding the kinetics. After this is accomplished 
and using the information about the kinetic rates, the model will be modified 
to predict data for other experimental systems. This will further define the 
reaction kinetics of three iodides (i - CjFyl, t - C^Fgl, and n - C^Fgl) and 
allow the establishment of scalability of the solar-pumped iodine laser for 
space-power applications. 
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Table I- Reaction rate coefficients found by fitting the data shown in figure 
1, 2, and 3. The reaction rates found are within previously 
published values. 


Reaction rate coefficient, (cm ) /sec 


Symbol Reaction 


R = i-C 3 F 7 


R = n-C^Fg R - t-C^F g 


K 1 R + I* = RI 
K 2 R + I * RI 

k 3 r + r = r 2 

K 4 R + RI = R 2 + I 

K 5 I 2 + R = RI + I 


3.48 x 10’ 12 
1.38 x 10" 11 
4.15 x 10" 12 
7.66 x 10' 17 
6.27 x 10‘ 13 


7.75 x 10 
2.88 x 10 
6.05 x 10 

7.76 x 10 
1.94 x 10 


-14 

-11 

-12 

-17 

-13 


1.66 

1.49 

3.56 

7.58 

2.34 


x 10 
x 10 
x 10 
x 10 
x 10 


-12 

-11 

-14 

-17 

-12 


C 1 I* + I + RI = I 2 + RI 1.09 x 10" 33 

C 2 I + I + RI = I 2 + RI 2.42 x 10" 32 

C 3 I* + I + I 2 = I 2 + I 2 8 x 10" 32 
C 4 I + I + I 2 = I 2 + I 2 3.8 x 10" 30 


1.11 x 10‘ 33 

1.09 x 10' 33 

2.76 x 10" 32 

2.28 x 10 -32 

8 x IQ* 32 

8 x 10 -32 

3.8 x 10' 30 

3.8 x 10 -30 


Q 1 I* + RI = I + RI 1.08 x 10 

Q 2 I* + I 2 = I + I 2 4.96 x 10 


1.17 x 10" 15 8.45 x 10 

4.96 x 10" 11 4.96 x 10 


73 


List of Figures 



Figures 1, 2 and 3. Theoretical predictions compared with experimental data 

for a flashlamp pumped laser oscillator with cylindrical 
geometry . 

Figures 4 through 12. Theoretical predictions compared with experimental data 

for the power output for a flashlamp pumped laser 
oscillator with cylindrical geometry for three different 
iodides and three pressures. 

Figures 13 through 21 Using the kinetic rate coefficients given in Table I for 

different iodine gases. The derivative with respect to 
time of the inversion density versus time is plotted for 
the different reactions given in Table I. The 
derivatives are renormalized to the sum of the 
photodissociation rates for RI and Ig and the pulse 
power output is normalized to its peak. The peaks of 
the derivatives and the power are given in the key for 
the different reactions. 
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Figure 
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Figure 2 
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Figure 
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Figure 





GAS I-C3F7I 
REACTION PEAK VALUE 



Figure 13 


GflS I-C3F7I 
REACTION PEAK VALUE 



Figure 14 


GAS I-C3F7I 
REACTION PEAK VALLE 



Figure 15 




GfiS N-C4F9I 
REACTION PEAK VALUE 



Figure 16 




Figure 17 



GAS N-C4F9I 
REACTION PEAK VALUE 



Figure 18 
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Figure 19 



GAS T-C4F9I 
REACTION PEAK VALUE 



Figure 20 



GfiS T-C4F9I 
REACTION PEAK VALUE 



Figure 21 




